Stability of superlubric sliding on graphite 
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Recent AFM experiments have shown that the low-friction shding of incommensurate graphite 
flakes on graphite can be destroyed by torque-induced rotations. Here we theoretically investigate 
the stability of superlubric sliding against rotations of the flake. We find that the occurrence of 
superlubric motion critically depends on the physical parameters and on the experimental condi- 
tions: particular scan lines, thermal fluctuations and high loading forces can destroy the stability of 
superlubric orbits. We find that the optimal conditions to achieve superlubric sliding are given by 
large flakes, low temperature, and low loads, as well as scanning velocities higher than those used 
in AFM experiments. 
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I. INTRODUCTION 

Recent years have witnessed a surge of interest in un- 
derstanding the microscopic origin of friction as a result 
of the increased control in surface preparation, the de- 
velopments of local probes like the Atomic Force Mi- 
croscopes (AFM) and Scanning Tunneling Microscopes 
(STM) and due to the interest for possible applications 
in nanotechnology. One of the goals of this research is 
to understand whether extremely low friction can be ob- 
tained by an appropriate choice of the sliding conditions. 
This paper examines theoretically the sliding of graphite 
flakes on a graphite substrate, one of the prototype sys- 
tems in this field. For this system, it has recently been 
shown that the low friction 'superlubric' sliding reported 
previously for flakes with incommensurate contact with 
the substrate [1 is always destroyed by rotations of the 
sliding flake [2 , leading to a locking in a commensurate 
state with high friction and slip-stick behavior. Numeri- 
cal simulations [2 carried out for the experimental condi- 
tions (extremely low velocities, about 30 nm/s) confirm 
this finding. It is intriguing to ascertain whether there 
might be conditions that avoid the rotation and locking 
in the high-friction commensurate orientation. 

Some important concepts of friction at the atomic scale 
are based on the Frenkel Kontorova (FK) model [84 that 
describes the sliding surface as a harmonic chain of lattice 
spacing a in interaction with a rigid periodic substrate 
with period h. For incommensurate values of the ratio 
a/6, Peyrard and Aubry [4] have shown that, below a 
critical value of the coupling to the periodic potential, the 
chain can be displaced on the substrate by an infinitesi- 
mally small force, namely the system displays a vanish- 
ing static friction force. Later, Shinjo and Hirano[5 pre- 
dicted that for incommensurate contacts also the kinetic 
friction would vanish and called this effect superlubricity. 
The experimental STM [6 and AFM studies [P showing 
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a drop of the friction force in going from commensurate 
to incommensurate contacts seemed to confirm the pre- 
diction of superlubricity. Theoretical work [7 has shown 
that the prediction of frictionless sliding also at high 
velocities of Ref. [6j is oversimplified and does not ap- 
ply in general, although dissipative mechanisms become 
less and less effective in the limit of vanishing velocities. 
Moreover the term superlubricity has been criticized in 
several papers [8^, ^ because it suggests a transition to 
zero friction which can be compared to superfluidity or 
superconductivity, whereas there is no threshold value 
of the velocity below which the kinetic friction vanishes. 
Nevertheless, the term superlubricity has become very 
popular and is used to describe low friction in the qua- 
sistatic limit accessible by AFM. 

Here we study the driven dynamics of a finite graphite 
flake on a graphite surface. The flake-surface interaction 
is modeled with a realistic static potential but vibrations 
of the flake are not taken into account and those of the 
substrate are represented by an effective friction coeffi- 
cient proportional to velocity. This defines a determin- 
istic non-linear dynamical system with four degrees of 
freedom that can be studied by numerical simulations 
and approximate analytical models, allowing us to study 
the stability of superlubric sliding. 

For a commensurate contact, we always find a stick- 
slip behavior with high friction. Conversely, for an in- 
commensurate contact we find two types of qualitatively 
different behavior. After an initial short period, the flake 
either rotates and locks into a commensurate orientation 
or it remains incommensurate and slides with extremely 
low friction. This behavior is critically dependent on the 
initial conditions, as expected for a strongly nonlinear 
problem. A simple dynamical system which captures 
the essential physics and for which the stability analysis 
can be done analytically explains the observed behavior. 
We then examine by numerical simulations the stability 
of the periodic orbits corresponding to incommensurate 
sliding against thermal fluctuations and other perturba- 
tions. 

In Sec. ini we describe the model of the structure and 
interactions and the details of the numerical simulations. 
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FIG. 1: Top view of the geometry of a graphite flake of 
24 atoms on the substrate, in a commensurate orientation 
(left, mismatch angle = 0) and incommensurate orientation 
(right, = 30°). The open circles represent substrate atoms, 
while the closed circles are flake atoms. The scan lines used 
in this paper are along the x-axis and shown from top to 
bottom: scan line 1 (solid line), 2 (dashed line), 3 (dotted 
line), and 4 (dot-dashed line). The scan lines are separated by 
a distance a/4. Due to the symmetry of the lattice, the range 
between scan lines 1 and 4 fully describes all scan lines in this 
direction. The scan line at a distance a/4 below scan line 4 
is again equivalent to scan line 3. Note that in a symmetric 
hexagonal flake, the center of mass does not correspond to 
the position of an atom. 



In Sec. Ill we show that periodic orbits corresponding to 
either commensurate or incommensurate sliding appear 
for different initial conditions. In Sec. |IV| we propose a 
simplified model for which we can perform analytically 
the stability analysis of these orbits. The robustness of 
the stability of periodic orbits against different types of 
perturbations is presented in Sectior|V| Finally we con- 
clude with a summary and perspectives in Section VI 



II. MODEL 

We study the dynamics of rigid graphite flakes, lying in 
the x—y plane parallel to the substrate as shown in Fig.[l] 
Atoms are kept at the equilibrium inter- atomic spacing 
a = 1.42 A in a hexagonal lattice for both the flake and 
substrate. By changing the orientation of the flake onto 
the hexagonal substrate the contact is either commen- 
surate (Fig. Ill left) or incommensurate (Fig. fll right). 
We consider only rotations around the z axis that keep 
the flake parallel to the substrate. The center of mass 
of the flake is pulled along the indicated scan lines by a 
support moving at constant velocity Vg = ('Us,0,0). The 
flake therefore has 4 degrees of freedom: the coordinates 
of the center of mass, r = (x, y^ z) and the orientation <p. 
The corresponding velocities are v = (vx^^y^^z) and uj. 
The phase space has 8 dimensions. 

We calculate the force and the torque acting on the 
center of mass from the interaction that each atom in 
the flake has with each atom in the substrate. The to- 



tal potential energy of the flake due to interactions with 
atoms of the substrate can be written as 



V{r,^)=Y,Y.^c{\r^-R, 



(1) 



where i goes over all flake atoms, and j over all sub- 
strate atoms and Vc{r) is the interaction between one 
flake atom and one substrate atom at distance r. The 
positions of the substrate atoms Uj = {Xj^Yj^ Zj) are 
given by a hexagonal lattice, and the positions of flake 
atoms Ti = (x^, yi^Zi) are functions of the position of the 
center of mass r = (x^y^z) and of the orientation angle 
(j) (see Fig. fl]). In the simulations described in this pa- 
per, we use the atom-atom interaction potential V^^(r) 
of Ref. [10 that describes non-bonded interactions of car- 
bon. The potential has a range of 6 A. 

The support representing the AFM cantilever drives 
the flake, with a force given by 



Fsir,t) 




-Fload 



(2) 



where t is the time, (xs,ys,^s) = (^s(O) + Vst,ys{0), z^) 
is the position of the support, c (= 1 nN/nm) is the 
coupling constant between the support and the center 
of mass of the flake, and Fioad is the load force in the 
negative z direction. The coupling to the phonon modes 
of the substrate can be modeled by a viscous friction term 
that dampens the motion of the flake, with a force and 
torque given by 



^f(v) = -7MV 



(3) 
(4) 



where M is the total mass of the flake, / is the moment of 
inertia for rotations around the center of mass along the 
z-axis, and 7 (= 1/ps) is the viscous friction constant. 
Note that for a rigid flake, the damping of the linear ve- 
locity directly determines the damping of both the centre 
of mass and the rotation. 

The equations of motion are: 



Mr. 
I4>. 



dV{r, 



dr 

dV{r, 



+ F,(r,i) + Ff(v), (5) 

+ TfH. (6) 



The rotational symmetry of the flake implies that 

V{r,<p) = V{r,^ + <P). (7) 

and the periodicity of the substrate gives 

V{r,^)=V{r + a,^), (8) 

where a is any vector which generates a translation under 
which the lattice is invariant. The flake-substrate system 
also has symmetry for reflections in the ?/z-plane 



V(r,^) = V{{-x,y,z),7r- 



(9) 



In our numerical simulations, we solve the equations of 
motion using the velocity- Verlet algorithm with damping 
and whenever the temperature is nonzero, a Langevin 
noise term is added [21 HI] • 



III. PERIODIC ORBITS 

The solutions of Eqs. ( 5| and ([g]) at T = are strongly 
dependent on the initial conditions, due to the nonlin- 
earities of the interaction forces. In Fig. [2] (top left) we 
show two trajectories obtained for exactly the same con- 
ditions (same load, support velocity, and scan line) apart 
from different initial angular velocity. We can see that 
starting from an orientation near the incommensurate 
orientation, cj) either drops to the commensurate (/) = 
value, or oscillates around approximately 26°. A similar 
trajectory on another scan line converges to 30°. The 
orientation converges to a stable value within a few lat- 
tice periods. This result shows that several periodic or- 
bits may be stable. The corresponding behavior of x(t), 
shown in Fig. |2] (top right), for the commensurate case 
^ = is step-like, which is typical of stick-slip motion. 
For the incommensurate cases (j) = 26°, 30°, the flake 
follows the support closely. The difference between com- 
mensurate and incommensurate orbits is also evident by 
looking at the trajectory in the x?/-plane, shown at the 
bottom left of Fig. [2J In the case of </> = the cen- 
tre of mass jumps quickly from one lattice site to an- 
other, where it performs some oscillations before jump- 
ing again. The incommensurate motion at the same scan 
line is smoother and the orbit at (/) = 30° performs a 
regular zig-zag motion. The lateral force, also displayed 
in Fig. [2] (right bottom), which shows stick-slip motion 
for the commensurate trajectory, drops for (j) = 26° and 
(/) = 30° to an average friction force close to that of a flat 
surface {jMvs = 0.0153 nN), 0.0278 nN and 0.0316 nN 
respectively. The friction of the commensurate flake, by 
comparison, is large, 0.1018 nN. 

In rare cases, particularly at very high load, where the 
nonlinearities are increased, periodic trajectories with a 
period longer than one lattice period as well as chaotic 
trajectories exist. Examples of a period 6 periodic orbit 
and a chaotic orbit are displayed in Fig. |3] Neverthe- 
less, even in these trajectories, the orientation remains 
roughly constant. 

In Fig. [4] the stable periodic orbits are plotted as a 
function of ?/s, ranging between scan line 1 and 4, for 
the system of Fig. [2] The commensurate periodic orbit 
at ^ = is always stable, regardless of the scan line. 
Between scan lines 3 and 4 the incommensurate orbit 
at ^ ~ 26° becomes unstable, and the one at (/) ~ 30° 
becomes stable. 

As the number of atoms increases, the interaction with 
the substrate becomes more complicated and the number 
of periodic orbits increases. For square flakes on a square 
lattice, the number of periodic orbits increases linearly 
with the diameter of the flake [12 . In Fig. pj bifurcation 
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FIG. 2: (Color online) Three typical trajectories for a 24- 
atom flake subjected to Fioad = 20 nN, Vs = 32 m/s. All three 
converge to stable periodic orbits at approximately constant 
(j) ^ (t)Q. The trajectories converging to 00 ~ 0° and 00 ~ 26° 
are for scan line 3, but have different initial angular velocity, 
and the trajectory converging to 00 ~ 30° is for scan line 4. 
From left to right and top to bottom: (a) the mismatch angle 
as a function of time, (b) the position as a function of time 
once the trajectories have converged to the periodic orbits, 

(c) the trajectories on the surface for the same interval, and 

(d) the friction force. 



(a) 



Mi'm Tlj ^1 Tlj Til Til 




(b) 



100 150 200 250 300 

tips) 




50 100 150 200 250 300 

tips) 



FIG. 3: Examples of more complicated trajectories of a 24- 
atom flake: (a) a periodic trajectory with a longer period, in 
this case 6 lattice periods, for Fioad = 30nN at scan line 3 
and (b) a chaotic trajectory for for Fioad = 40nN at scan line 
1. All other conditions are the same as for the trajectories 
plotted in Fig. |2] 



diagrams similar to Fig. [4] are shown for flakes of different 
sizes. The number of stable periodic orbits increases. Ad- 
ditionally, there is a switch-over region around scan line 
3, where the stable incommensurate orbits become un- 
stable, and the unstable incommensurate orbits become 
stable. 

Experimentally [2] it was reported that the superlubric 
behavior of flakes of approximately 100 atoms lasted for 
about 40 scan lines or a distance y^ of about 7A, about 
5 times the distance between scan lines 1 and 4. This 
compares very well with the results for A^ = 96, where 
we see that starting, for instance on scan line 4 and mov- 
ing towards scan line 1, the flake rotates from the sta- 
ble periodic orbit at 30° to the one at 23°. After this. 
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FIG. 4: A bifurcation diagram of the stable periodic orbits as 
a function of the parameter y for N — 24: and Fioad = 20 nN. 
The data was obtained by doing a large number of simulations 
with a wide range of initial conditions. The plotted points are 
the set of final angles. Clearly visible between scan lines 3 and 
4 are the points at which the 00 ~ 26° periodic orbit becomes 
unstable and the 00 ~ 30° periodic orbit becomes stable. 



due to the symmetry of the lattice, the scan moves back 
from scan line 1 to scan line 4, decreasing the mismatch 
angle to 19° (which has lower energy than 30°). After 
3a ~ 4.3 A distance in the y direction, the flake locks 
in the commensurate = state. In absence of ther- 
mal fluctuations the decay to the commensurate state is 
a geometric effect, depending only on the structure of 
the interaction. Each periodic orbit leads to a different 
friction force, and so the observation of steps in the fric- 
tion force going from one scan line to another, could be 
related to the size and symmetry of the flake. 
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FIG. 5: The plot of Fig. HI repeated for (a) various 6-fold 
symmetric flakes of (b) 54, (c) 96, (d) 150, and (e) 216 atoms. 
Larger flakes have more stable periodic orbits. The lone point 
in the bifurcation diagram for N — 9Q near scan line 4 at 
00 ~ 12° indicates that the 0o ~ 12° periodic orbit is still 
stable there, but has such a small basin of attraction that the 
spacing between the initial conditions used to calculate this 
bifurcation diagram is not fine enough to detect it. 



this section. 



IV. STABILITY ANALYSIS AND SIMPLIFIED 
MODELS 

Although we consider the flake as a rigid object with 
only four degrees of freedom, the system is still too com- 
plicated to perform the stability analysis analytically. 
However, a possible simplification is suggested by the 
shape of the potential energy ]/(r, 0). In Fig.[6J we show 
y as a function of x and for constant values of z given 
by the average of the values found in simulations with 
a load of 20 nN, and y defined by the four scan lines. 
One can see that the potential is a periodic function of 
X with an amplitude that depends on 0. Therefore, a 
good description of the system is provided by a simplified 
one-dimensional model with only two degrees of freedom: 
the position along the scan line, x, and the orientation, 
0. This model is fully described by the viscous friction 
coefficient 7, support velocity v^^ initial support position 
Xg , mass M, moment of inertia /, and a simplified poten- 
tial y(x,0). The essential dynamics of the system, the 
existence of commensurate and incommensurate sliding 
is preserved in the simplified model which we present in 



A. equations of motion 

We write the equations of motion of the simplified sys- 
tem as a dynamical system of first-order differential equa- 
tions. 



(10) 

MVr. = - '"' r''^^ - C{x - tVs - Xs 0) - jMVa; , (11) 
</. = W, (12) 

^- = -^^^-7/-. (13) 



dV{x, 4>) 
dx 



Moreover, the symmetries of V{v, (p) in Eqs. W^ im- 
ply that 



V{x,<j)) = V{x,- + ~ 
V{x,4>) = V{x + l,^ 

V{x,(t))=V{-X,Tl-^ 



(14) 

(15) 
(16) 



where / = aV^- 
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FIG. 6: (Color online) The potential energy V(r, 0) of a 24- 
atom flake as a function of the mismatch angle (f) and position 
X along the trajectory of the support, for constant z at the av- 
erage value belonging to a load of 20 nN, and y corresponding 
to scan lines (a) 1, (b) 2, (c) 3, and (d) 4. Due to the sym- 
metries of the system given in Eqs. (14} 16), the dependence 



on (/) is determined by the behavior between and 30° 



1. Specific potential 

A good representation of F(x, (j)) for a given scan line 
{y constant) is given by 



/27rr\ 
V{x, 4>) = U{(f) + W{4>) cos ( — J 



(17) 



where t/(^) and W{(j)) are both smooth functions that 
represent the average value of the potential energy and 
the amplitude respectively. 



The symmetries of the dynamics in Eqs. ( 14 - 16 ) imply 
that 



[/(0) = f/(- 

w{4>) = w{- 






(18) 
(19) 



In turn, these equations imply that U and W have ex- 
trema in = 0o = 0, 7r/6. In Figs. [tI and [sj U and W, 
are shown for flakes of 24 and 216 atoms. It is evident 
in Fig. [7| that there is an extremum of both U and W 
at (/) = 0. The structure of the other extremum, close to 
30° can be seen from the four enlargements. Besides the 
extremum at ^ = 30° for all scan lines there is another 
extremum of both U and W at about 26°. In Fig.jS) for a 
larger flake, U and W have more extrema, but they still 
coincide. In Ref. [12, it is shown for a simpler system, 
square flakes on a square lattice, that this is a general 
property: for square flakes on square lattices the extrema 
of U and W at any orientation coincide approximately for 
all flake sizes. 



Since the torque, given by Eq. (13), vanishes for the 
values of (j) that give extrema of U and W and a; = 0, 
these conditions define a two-dimensional invariant man- 
ifold of the dynamics. The number of extrema of U and 
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FIG. 7: The (a) offset [7(0) and amplitude (b) 1^(0) of 
the potential 1/ (x, 0) as a function of (j) for the same case 
displayed in Fig. [6] The region near = 30° is enlarged 
separately for scan lines (c) 1, (d) 2, (e) 3, (f) 4. U and W 
were obtained from a Fourier transform of V with respect 
to X over 492 points for each 0. The extrema of U and W 
coincide at = 00 , which implies the existence of an invariant 
manifold = (/)o,cl; = with 00 = 0°, 26° or 30°. 



W and consequently the number of invariant manifolds 
grows with the size of the flake. 



B. stability 

We consider a general potential V(x, 0) which has an 
invariant manifold at (j) = 00, i.e. 



SV(x,0) 



= 



(20) 



for all X. Now that we have identified the invariant man- 
ifold, we consider the dynamics in its vicinity in order to 
study the stability. 

Near the invariant manifold, the torque is small, and 
so the time scales of and cj, (Eqs. ([l2| and (13)) are 
much longer than those of x and Vx (Eqs. (10) and (11)). 
Because of this, for the purpose of investigating the sta- 
bility of the dynamics near the invariant manifold, the 
torque can be replaced by its time average. Note that 
this separation of time scales is only valid near the in- 
variant manifold, namely if remains close to 0o and uj 
is close to 0. 
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FIG. 8: The (a) offset [/(0) and (b) amplitude 1^(0) for a 
flake of 216 atoms. The extrema of U coincide with the max- 
ima of VK, and the nodes of W correspond to a constant value 
of U. There are more extrema than for the flake of 24 atoms, 
and therefore more stable and unstable incommensurate pe- 
riodic orbits. 



If the manifold is stable, then initial conditions close to 
it converge towards it. We therefore consider the growth 
rates of small perturbations 5(j) and Suj of (j) and cj, the 
Lyapunov exponents. From Eqs. (12) and (13) we find 



S(j) = Suj ^ 
I Suj = —Sd) 



d /dV{x,i 



(21) 
■ jISuj . (22) 



The time average can be interchanged with the deriva- 
tive with respect to (j) because perturbations in x and (j) 
decouple to first order. One may write 



Suj 







1 



A 



1 / d'^Vix,< 



Suj 



Suj 



(23) 

(24) 



As the matrix A is constant, the Lyapunov exponents 
associated with perturbations in (j) and uj are simply equal 
to its eigenvalues. 



--7±- 
2^ 2\ 



4 I dW{x,(j)) 

1\ w 



(25) 



0=00 I 



The invariant manifold is stable if all (in this case 2) 
Lyapunov exponents associated with perturbations of it 
have real components smaller than 0. 



As the real component of the square root in Eq. (25) is 
positive or 0, A_ < A+ is the smallest Lyapunov exponent 
(i.e. has the smallest real component). For stability anal- 
ysis it therefore suffices to consider A+ . If the argument 
of the square root in Eq. (25 ) is smaller than 7^, then the 
real components of both A_ and A+ are negative. This 
is the case if 



o>2V(x, 



>0 



(26) 



i.e., the time-average of the potential energy must be at 
a minumum. 



Using Eq. (17), Eq. (26) can be rewritten to read 

d^U{(l)) d^W{(l)) 



f2TTx\ 



cosl^l ) >0. 

(27) 



The stability thus depends on U and W, and how much 
time the particle spends near the minima of the potential, 
where the cosine is negative. 

In stick-slip motion, the particle spends most of its 
time in the minima of the potential, i.e. where the cosine 
is smaller than (see Fig. |2|. If the motion is truly 
superlubric, then the particle spends about the same time 
in the minima as it does in the maxima. If the motion is 
nearly superlubric, then the particle spends most of its 
time in the minima. Hence, for realistic cases, (cos)t < 0. 

If the offset of the potential, t/(^), has a minimum at 
(j)Q it contributes positively towards the stability. Sim- 
ilarly, if the amplitude W{(j)) is at a maximum at ^o, 
because the first derivative is multiplied by a negative 
number, (cos)^, it enhances the stability. A minimum of 
U and maximum of W therefore always lead to stabil- 
ity, whereas a maximum of U and minimum of W always 
leads to instability. If both are at a maximum, or both 
are at a minimum at ^o^ then the stability is not directly 
obvious. 



comparison with simulations 



= 00 I 



The analysis of Sec. IV B compares very well with the 
results of numerical simulations at T = K. The stabil- 
ity of the commensurate and incommensurate states can 
be determined by looking at the behavior of the average 
potential energy U and amplitude W, shown in Figs. [7| 
and El 

We examine first the 24-atom system of Figs. [2] and [4j 
for which U and W are reported in Fig. [7| For scan line 1 
and 2, the minimum of /7 at (/) = 26° coincides with a 
maximum of W ^ and is therefore stable. This is consis- 
tent with the simulation results for scan lines 1 and 2, 
shown in Fig. ^ where we see a stable orbit at 26° . At 
these scan lines, for (j) = 30° U has a maximum and W 
has a minimum, leading to instability. At ^ = 0, there 
is a maximum in [/, but also in W . However, the second 
derivatives of U and W are very nearly the same apart 
from the sign, and (cos)^ increases with decreasing ampli- 
tude, so {d'^V/d(j)^)t is positive, and the incommensurate 
state is stable. 

For scan lines 3 and 4, at 0° the minimum of U co- 
incides with a maximum in W ^ leading to a stable com- 
mensurate state. Similarly, at scan line 4, the incom- 
mensurate state at ^ = 30° is stable, while the state at 
(j) — 26° is unstable. For scan line 3, the stability of the 
incommensurate states is more complicated, as U and 
W both have maxima around 6 = 26° and minima at 



(I) = 30°. However, the second derivatives of U and W 
for both states are approximately the same, with oppo- 
site sign. Additionally, the amplitude for both states is 
approximately the same, so (cos)t should be the same 
as well. The crucial quantity for stability, {d'^V/d(j)^)t^ 
should therefore be nearly the same for the two states, 
except for the sign, which is opposite. One of the incom- 
mensurate states is therefore stable, while the other is 
unstable, though from U and W it is not directly clear 
which is which. In Fig. [2J the stable incommensurate 
state for scan line 3 is shown at = 26° and for scan 
line 4 at = 30°. The existence of a switch-over scan 
line can be seen in the simulation results in Fig. [4j and is 
clearly critical for all sizes, as shown in Fig. [5] Its exis- 
tence for any flake size can be demonstrated analytically 
for square flakes on square lattices [12 . 

In Fig. [8J U and W are plotted for a larger flake of 
216 atoms. Because of the larger size of the flake, U 
and W have more extrema and therefore there are more 
periodic orbits. The stable periodic orbits in the simula- 
tions, shown in Fig. [s] (bottom right), coincide with the 
extrema. Their stability is also consistent with calcula- 
tions based on U and W. 



V. ROBUSTNESS OF THE SUPERLUBRIC 
SLIDING 



The analysis presented in Sec IV shows that incom- 
mensurate (superlubric) sliding may exist. However, the 
existence of stable incommensurate periodic orbits does 
not necessarily mean that they can be easily observed in 
experiments. The conditions which lead to stability may 
not be experimentally accessible. Furthermore, the sta- 
bility may be very weak, causing very slow convergence 
towards the periodic orbit, or the basins of attraction 
of the incommensurate periodic orbits (the set of initial 
conditions that converge towards them) may be small. 
In this section we examine separately the robustness of 
the incommensurate superlubric solutions against several 
types of perturbations. 



A. temperature 

In Fig. [9J we show slices of the phase space which con- 
tain the stable incommensurate periodic orbits of the full 
three-dimensional system. For each scan line, we inves- 
tigate the basin of attraction by performing numerical 
simulations at T = and looking at the asymptotic state 
of the flake as a function of the initial orientation and 
angular momentum. If the basin of attraction is small, 
the periodic orbits can easily be destroyed by thermal 
fluctuations, which bring the system outside the basin of 
attraction of the incommensurate orbits, and into that 
of the commensurate orbit. The range of initial angu- 
lar velocities plotted in Fig. Mis 3y^k^T^jM^ where T^ 
is room temperature, 293 K, and /cb is Boltzmann's con- 



stant. This is roughly the range that is thermally ac- 
cessible at room temperature. The basin of attraction 
of the incommensurate periodic orbits is smaller than 
this range, indicating that at room temperature thermal 
fluctuations may perturb the incommensurate state suf- 
ficiently to cause it to decay to the commensurate state, 
which has lower energy. Especially scan line 3, with its 
weak stability and scan line 4, at which the incommensu- 
rate state only has a small basin of attraction (as is shown 
in Fig. [9]), are very sensitive to thermal fluctuations. 

To examine the effect of temperature explictly we con- 
duct Langevin simulations with temperatures ranging 
from 5 to 300 K. Starting from initial conditions on the 
incommensurate periodic orbit, simulated systems were 
subjected to thermal fluctuations for a period of about 
100 lattice periods and the final angle was recorded. The 



results are plotted in Fig. 10 At scan line 3 the incom- 
mensurate state is the least robust against temperature 
and the incommensurate state decays already at 5 K. 
However, thermal fluctuations are not the only source of 
energy in this system, because the moving support drives 
the flake at velocities that are not negligeable compared 
to thermal velocities, therefore supplying amounts of en- 
ergy significant compared to k^T. The effect of temper- 
ature is thus possibly overestimated in these simulations. 



B. scan line 

From the size of the basins of attraction in Fig. [9] and 
the robustness against thermal fluctuations, displayed in 
Fig. [lOJ it can be seen that the robustness of the incom- 
mensurate states in this system depends strongly on the 
scan line. For the system in the figures, the incommen- 
surate periodic orbit is the least robust for scan lines 3 
and 4. From Fig. [7| it can be seen that the minimum of 
U near (p = 30° is shallow and the amplitude W, espe- 
cially in the case of scan line 4, is small. The latter is a 
consequence of the symmetries of the hexagonal lattice. 

As discussed in Sec. |III[ the different stability and in- 
stability of the incommensurate states at different scan 
lines can lead to the disappearance of superlubricity after 
an initial superlubric period in experiments which explore 
more than one scan line. Additionally, the weak stability 
of the incommensurate states, and associated low robust- 
ness against thermal fluctuations, near scan lines 3 and 4 
makes superlubric states less likely to persist in such ex- 
periments. 



C. flake size 

As the number of atoms in the flake increases the mo- 
ment of inertia increases with N'^. This means that the 
orientation and angular velocity of larger flakes are less 
sensitive to thermal fluctuations and other disruptions. 



By comparing Fig. [TT] to Fig. [TO] we see that the incom- 
mensurate periodic orbit of the flake with 216 atoms is 







oatj = n{Vpp) 



inmmmenGuretB 



tttttt 

-12D -SD -BD -an 



final 1^ (dsgrBBs) 



eammenGurete 



FIG. 9: (Color online) Cross sections of the phase space, including the basin of attraction of the incommensurate stable 
periodic orbits, which are at ^ 26°, cj = 0, for scan lines (a) 1, (b) 2, (c) 3, and ^ 30°, cj = for (d) scan line 4. The flake 
has 24 atoms and Fioad = 20 nN, Vs = 32 m/s. The final state of the flake is plotted as a function of the initial orientation 
and angular momentum. The initial position and velocity have been chosen such that the stable incommensurate periodic 
orbit intersects with the cross section, in cj = 0. The colours indicate to which periodic orbit the initial conditions converge: 
red incommensurate 00 G (0°,30°], blue incommensurate (j)o G [30°, 60°), purple with blue commensurate 00 = 60°, purple 
incommensurate 00 G (60°, 90°], black incommensurate 0o G [90°, 120°), red with black commensurate 0o = 120°, green with 
red commensurate 0o = 0°, green incommensurate 0o G [— 30°,0°), cyan incommensurate 0o G (—60°,— 30°], yellow with 
cyan commensurate 0o = —60°. yellow incommensurate 0o G [—90°, 60°). grey incommensurate 0o G ( — 120°, —90°]. The 
incommensurate periodic orbit at 0o ~ 30° is indicated in blue, as it visits both (0°,30°] and [30°, 60°) in one period. 



more robust against thermal fluctuations, and survives 
even at room temperature. 

It is interesting to note that Bonelli et al.[13 , who 
consider flexible graphite flakes, found that larger flakes 
interact more weakly with the substrate than one would 
expect from rigid flakes. At the edges, the flake bends 
towards the substrate, and thus the atoms at the edge of 
the flake dominate the interaction. However, in Ref. [13], 
no analysis of the stability of superlubricity was possi- 
ble, as the coupling between the cantilever and flake was 
chosen in such a way as to impose a preferred orientation. 



D. support velocity 

At high support velocity, the motion of the flake is 
less sensitive to the detailed structure of the substrate. 
The dynamics in the y and z direction are relatively fast 
compared to the dynamics of the rotation, and therefore 
their effects on the orientation of the flake average out. 
At lower support velocities, motion in the y and z di- 
rection becomes more relevant and can reduce the size 
of the basin of attraction of the incommensurate peri- 
odic orbits, or even destroy the stability completely. In 
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FIG. 10: The final orientation of a 24- atom flake (mapped 
onto the interval [0°,30°]) which was initially in the stable 
incommensurate periodic orbit is plotted as a function of tem- 
perature after a long, but finite time, with Fioad = 20 nN, Vs = 
32 m/s for scan lines (a) 1, (b) 2, (c) 3, and (d) 4. For every 
temperature, 250 realisations are plotted. Enough time has 
elapsed for the system to decay to the static state distribution. 
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FIG. 12: (Color online) The (a) orientation 0, positions (b) x 
and (c) y, and (d) friction Fs as a function of support position 
Xs for various support velocities and A/" = 24, Fioad = 20 nN, 
scan line 2. As the velocity decreases, the fluctuations in 
and y increase and the system behaves less one-dimensionally. 
For sufficiently low Vs, the system can no longer be described 
by the simplified model. 




(a) 



[ [□□□□□□□□□□□□□p 



10 - 
> 



1 2 3 

y^ (scanline) 



(b) 



I aB H uu n uu Bii 



iHQDQHCiQEiQCiniDnHan^^ 



2 3 

y^ (scanline) 



FIG. 11: The plot of Fig. [To] for scan line 2 repeated for 
flakes of (a) 96 and (b) 216 atoms. The stable periodic orbits 
of large flakes are more robust against temperature, because 
the moment of inertia grows as N^ . 



Fig. 12 trajectories are plotted for the same flake at dif- 
ferent support velocities. As the velocity decreases, the 
flake becomes more sensitive to fluctuations and there- 
fore (j) (top left) and y (bottom left) fluctuate more. At 
sufficiently low velocities, the incommensurate periodic 
orbit is no longer stable, and the flake rotates to the com- 
mensurate orientation with stick-slip motion (top right) 
and high friction (bottom right). A stronger coupling 
between the flake and cantilever would reduce the fluc- 
tuations in the y direction, and allow the stability of the 
incommensurate state to persist to lower support veloci- 
ties. 



E. load 

The load force exerted by the cantilever on the flake 
pushes it into the substrate. This affects not only the 
corrugation, but also the shape of the potential to which 




1 2 3 

y^ (scanline) 



2 3 4 

y^ (scanline) 



FIG. 13: The plot of Fig. |4] repeated for load force equal to 

(a) nN, (b) 10 nN, (c) 30 nN, and (d) 40 nN. 



the flake is subjected. Consequently, for different load 
forces, the behavior of U and W is different, and so the 
stability of incommensurate periodic orbits may change. 
In Fig.[T3]bifurcation diagrams similar to the one in Fig. [4] 
are shown for different load forces. When the load is very 
high, the interaction between the flake and substrate is 
changed qualitatively, and for the region near scan lines 3 
and 4 the incommensurate periodic orbit disappears. The 
simulations of Ref. [13] were performed using load forces 
of about 100 nN, and indeed, no superlubric behavior was 
observed. Bonelli et al. also performed a few simulations 
at lower loads for 24 atom flakes, but imposed mismatch 
angles near 0° and 15° on the flake, thus eradicating the 
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FIG. 14: (Color online) The cross section for scan line 3 in 
Fig. [9] repeated with the commonly used 2D interaction po- 
tential. The basin of attraction is different in shape and size. 



tential used here, as can be seen from Fig. [13 



However, in a fully three-dimensional problem, the ef- 
fect of load is not simply a rescaling of the amplitude. We 
compare our results obtained with a three-dimensional 
potential to the one obtained with the two-dimensional 
potential of reference [2 . We find (Fig. Il4| that the cross 
section has qualitatively the same features, but a signifi- 
cantly different size of the basin of attraction. 



G. symmetry of the flakes 



In experimental conditions, it cannot be guaranteed 



that the flakes are exactly hexagonal. In Fig. 15 the bi- 
furcation diagrams of Figs. |4] and [5] have been repeated for 
three- fold symmetric flakes of two different sizes. The re- 
sults are similar to those of the hexagonal flakes, though 
somewhat distorted. 



VI. CONCLUSIONS 
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FIG. 15: The plot of Fig. [I] repeated for (a) various 3-fold 
symmetric flakes of (b) 33 and (c) 69 atoms. 



incommensurate periodic orbits near 30°. 

At high load force periodic trajectories with periods 
longer than one lattice period and chaotic trajectories 
exist. Two such trajectories are shown in Fig. |3] These 
trajectories still have roughly constant orientation, be- 
cause the invariant manifold is still stable. It is the mo- 
tion on the invariant manifold itself that has a longer 
period or is chaotic. 



F. choice of potential 

Very often, for friction, the potential corrugation is 
represented as a two-dimensional profile in the xy plane. 
In this representation, the load can only be included by 
scaling the potential. Figs. [7| and [S] would therefore look 
the same but only scaled, regardless of load, which im- 
plies that the stable incommensurate orbits would remain 
stable for any load. This is not the case for the 3-d po- 



In this paper, we have examined the possibility of real- 
ising conditions for superlubric sliding without rotation 
and locking of graphite flakes on graphite. By means of 
a simplified analytical model, validated by our numerical 
simulations, we have shown that incommensurate peri- 
odic orbits with low friction can be stable. Furthermore, 
we have investigated the robustness of the superlubric 
sliding against changes in several conditions and quanti- 
ties: temperature, scan line, fiake size, support velocity, 
load, and asymmetry. 

Our results show that some scan lines, where the center 
of mass moves along a row of atoms of the substrate, are 
detrimental to the stability of superlubric sliding and lead 
to rotation of the flake as found in Ref. [2^ . Conversely, 
superlubric sliding is favored by larger flakes, higher ve- 
locities than in AFM, and low temperature. Our calcula- 
tions suggest that in an experiment where different scan 
lines are explored successively the locking would occur 
gradually via intermediate periodic orbits. For a flake 
of about 100 atoms, this should occur in 4 steps. As 
the friction force for each periodic orbit is different, this 
could perhaps be used as a method for characterising the 
flake. 
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